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Abstract 

We continue the semiclassical analysis, started in a previous paper, of the intertwiner sector 
of the flipped vertex spinfoam model. We use independently both a semi-analytical and a purely 
numerical approach, finding the correct behavior of wave packet propagation and physical ex- 
pectation values. In the end, we show preliminary results about correlation functions. 

PACS: 04.60.Pp 



1 Introduction 

In m [2] [3] jl] a new spinfoam vertex, now called the EPR or ELPR7 vertex, was introduced. This 
model, as well as other new models [5] [6] [7] , was born essentially in order to correct two problems 
of the Barret-Crane vertex |5]: the wrong intertwiner (non) dependence and the mismatching with 
the Loop Quantum Gravity [9 [10 [llj boundary states. As shown in [18], the BC model has a 
good semiclassical behavior in the spin sector, but not in the intertwiner sector; actually in the 
BC theory the intertwiner degrees of freedom are dynamically frozen. In particular, the diagonal 
components of the graviton propagator have the right scaling with distance [12J [13] [Pi] |15j [16] |17] . 
but the nondiagonal (the ones depending on the intertwiners) have the wrong scaling |18]|19j. 
An important open issue is to check if the new spinfoam models correct those pathologies. A 
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new technique to test the spinfoam dynamics, which is complementary to the calculation of n-point 
functions, was introduced in [50], and a partial answer was given. This technique is the propagation 
of semi classical wavepackets: as in ordinary quantum mechanics, if the theory has the correct 
semiclassical limit, then semiclassical wavepackets must follow the trajectory predicted by classical 
equations of motion. In [20 , the wavepacket propagation in the intertwiner sector was studied 
numerically, finding a surprisingly good semiclassical behavior. In brief, the authors considered the 
solution of discretized Einstein equations given by a single flat 4-simplex with boundary constituted 
by five regular tetrahedra. In the dual LQG picture this boundary is represented by a pentagonal 4- 
valent spinnetwork, labeled by ten spins and five intertwiners; but in order to have a semiclassical 
state one has to construct some (infinite) linear combination of spinnetworks of this kind. It is 
known from [2T] that basis 4-valent intertwiners with some choice of pairing can be superposed 
with gaussian weight to be able to catch the classical geometry: since in a quantum tetrahedron 
the angles do not commute, one has to consider semiclassical superpositions to peak all angles 
on the same classical value. In [20 the authors chose an initial state formed by four coherent 
intertwiners at four nodes, and they made the drastic approximation of taking the ten spins fixed 
to be equal to some jo. Then they calculated numerically its evolution, here called 4-to-l-evolution, 
that is its contraction with the propagation kernel of the flipped vertex spinfoam model. Classical 
Einstein equations impose the final state to be a coherent intertwiner with the same geometrical 
properties (mean and phase). They found good indications but, due to the very low jo's, they 
couldn't conclude much about the analytical properties of the evolved state. Here we conjecture 
the general behavior of the evolution at high jo's which is very well supported numerically. In 
fact, as we shall argue, the propagation is perfectly "rigid": four gaussian wavepackets evolve into 
one gaussian wavepacket with the same parameters, except for a flip in the phase. The phase of 
the evolved phase, and in particular its flipping, will have a major role when considering physical 
expectation values. The support to our conjecture is made in two independent ways: the first 
is semi- analytical and is based on a numerical result on the 15j-symbol viewed as a propagation 
kernel, and the asymptotic properties of the fusion coefficients already studied in [22] ; the second 
is purely numerical. The first has the advantage of giving a nice picture of the dynamics in terms 
of wavepackets evolving separately in the left and right £0(4) sectors, and it also pave the way 
for the completely analytical approach (to do this, one should have an asymptotic formula for the 
15j-symbol). We also explore the possibility of propagating three coherent intertwiners into two 
(we will refer to as the 3-to-2-evolution) , finding similar results. Then we present the results from 
another point of view, namely as intertwiner physical expectation values, finding that these are 
asymptotically the predicted ones. Though we use the drastic approximation of fixing all spins, 
we regard our results as a strong indication that the EPR model has the good semiclassical limit. 
In the end we present the numerical calculation of the intertwiner correlation function, finding a 
scaling law which is not the Newtonian one; we believe that this is due to the drastic approximation 
of freezing the spin variables in the boundary state, and not to pathologies of the model. 
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2 The flipped vertex fusion coefficients 



One of the key ingredients of the EPR dynamics is the branching function or fusion coefficients 
fl i , which are the coefficients of a linear map from SU(2) to SO(A) intertwiner vector space 
appearing in the definition of the vertex. In the case of 7 = 0, called flipped vertex, the fusion 
coefficients are defined by the evaluatiorP] of the spinnetwork in Fig. [I] They are symmetric under 




Figure 1: the flipped vertex branching function 

exchanging of ii and ir and non-vanishing only if the relation — ir\ < i < %l + Zr holds. In (52] 
the fusion coefficients were studied and an analytic asymptotic formula was given; thanks to this 
formula, we found the asymptotical action of the fusion coefficients on a semiclassical intertwiner. 
We resume briefly some properties. Consider a semiclassical regular tetrahedron, whose faces have 
(large) area jo; in LQG the tetrahedron, in particular its internal angles, is described by an SO(3) 
semiclassical intertwiner, defined as 

\k)sc = X^(0l*)> C 1 ) 

i 

where 

i>(i) = C(j ) exp (-^(i - i ) 2 + tfi) , (2) 

\i) is a basis intertwiner between four SO(3) irreps labeled by jo, C(jo) is a normalization constant 
and iq = 2jo/\ // 3 (i labels the virtual spin in some pairing). For simplicity, we will refer to ([T|, 
or to its components as "semiclassical intertwiner". We denote with t the imaginary unit, in 
order to avoid confusion with the intertwiner labels. The action of /| i (viewed as a map between 
intertwiner spaces) on a semiclassical intertwiner is given by 

g(iL, %r) = dim ( i )^,i ii V'(«)- (3) 

i 

1 The evaluation of a 3-valent spinnetwork is defined as the contraction (sum over magnetic numbers) of 3j'-symbols, 
one for each node, according to the pattern given by the graph. 
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We showed that, for large jo's 

g(i L ,i R ) * Oexp (- 4(i L - |) 2 - ^{i R - l f) 2 + L%(iL + iR)), (4) 

where C is an irrelevant normalization factor not depending on and i R at leading order in 1 /jo 
powers. Hence, asymptotically, the function g factorizes into left and right parts; we indicate 
them, with abuse of notation, g(i£) and g(i R ). The values of g(iL,i R ) are the components of an 
SO (4) ~ SU{2) x SU{2) intertwiner in the basis ir), which in the next we will call 50(4) 
semiclassical intertwiner. Also the converse holds: the asymptotical action of the fusion coefficients 
on an 50(4) semiclassical intertwiner is an 50(3) semiclassical intertwiner, i.e. 

dim(i L )dim(i i? ) f ih iR g(i L , i R ) ~ if>(i). (5) 



3 Propagation of intertwiner wavepackets 



The property Q gives a new picture of the dynamics in the semiclassical regime. Consider a 
boundary LQG state supported on the dual boundary of a 4-simplex (this state is labeled by ten 
spins and five intertwiners), and consider the simple case with all spins fixed to be equal to some 
jo, while the intertwiners are the semiclassical ones, namely: 
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*({jn},{i m }) oc J] 5 jndo Yl V(i 



(6) 



n=l 



where ift is the Gaussian Q. This state is the limiting case of a state peaked on the intrinsic and 
extrinsic geometry of the boundary of a regular classical 4-simplex as it becomes sharp in the spin 
variables [20 . We expect that, if four semiclassical intertwiners are given as initial state, then their 
evolution under the propagation kernel is again a semiclassical intertwiner (for an introduction and 
motivations to this idea see [20]), the state ([6j being peaked on a boundary solution of Einstein 
equations. Actually this propagation property should be true for a more correct boundary state in 
which the approximation of freezing the spins is not taken; nevertheless one can study the evolution 
and see if a positive result is obtained. If this is the case, we are strongly encouraged to believe that 
the same property holds in the general case. Coming back to details, define the 4-to-l-evolution as 



1*5 = 



dim(ii) . . . dim(u) W(i\ ... is) i/>(ii) ■ ■ ■ 



(7) 



Here 



w(h...i 5 ) = 



{inL}{inFl} 



Y[ dim(i nL ) dim(i nR ) 



.71=1 



15j(i 1L , i 5L ) 15j(i m , . . . , i 5R ) /? 



(8) 



n=l 
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Figure 2: On the left: modulus of the evolved state for the 4-to-l propagation performed by one 
15j (jo = 30). On the right: its real and imaginary (dashed) part. 



is the EPR flipped vertex amplitude evaluated in the homogeneous spin configuration (the ten spins 
are all equal to some fixed jo), and we have omitted the dependence of the 15j-symbol on the spins. 
The factor 4 . . . f^ L of Q is contracted in Q with four initial packets (making the sum 
over %\ . . . By Q, for large jo's this contraction gives four 5*0(4) semiclassical intertwiners, so 
the evolved state (|7| becomes 



^2 l5j(i 1L , i 5 L)g{iiL) • • • q{ul) 



^2 15 X i i«> • • • ' hR)g(hR) ■ ■ ■ g{uR) 



(9) 



We can see in the last expression the action of two 15j's separately on the left and right part (the 
expressions in square brackets). Those actions are interpreted as independent 4-to-l-evolutions in 
the left and right sectors, namely the evolution of the left an right part of four SO(4) semiclassical 
intertwiners, where the dynamical vertex is the 15j-symbol. By numerical investigations (Fig. [2j, 
it turns out that the final state of the right (left) partial evolution is the right (left) part of an 
SO(4) semiclassical intertwiner, with the phase flipped as compared to the incoming packets. For 
example, for the left part: 

0l(«5l) = ^ dim(ii L ) . . . dim(i 4L ) 15j(i 1L ■ ■ -hi) 9{hl) ■ ■ -9l{ul) — fife)- (10) 
Then ^ becomes 



ik) - ^2 dim i5L dim i5R 9{kL) g{kii) /, 

i5L,isR 



«5 



(11) 
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The last expression is the contraction between the fusion coefficients and an SO(4) semiclassical 
intertwiner. By ([H]), this gives an 50(3) semiclassical intertwiner: 

4>{^)~M^)- (12) 

While in [20 we expected only a conservation of mean values, and possibly a spread of wave packets, 
the precedent argument shows that the gaussian shape is conserved, together with its mean value 
and width, while the phase is flipped. The 3-to-2-evolution is defined similarly to the 4-to-l case, as 
the contraction between the flipped vertex and three initial semiclassical intertwiners. Numerical 
results about this type of evolution are discussed in section |5j 



4 Physical expectation values 

In this section we want to present the precedent results from another perspective, as results about 
expectation values of observables. By construction, the boundary state ([6]) is peaked kinematically 
on a semiclassical geometry. This should be also true in a dynamical sense, as it is peaked on a 
solution of Einstein equations. So consider the physical expectation value of an intertwiner on this 
boundary state: 

U '~ Z {j}>{i} W({j},{i})*({j},{i}) ' 1 1 

We expect this quantity to be equal to io for large jo's, if the dynamics has the correct semiclassical 
limit. Analogously, we can consider the expectation value of two intertwiners: 

,. _ J2 {3U > } W({j},W)m 2 n{j}A*}) . 

{l2) - E {j} , {i} wm,{i}M{j},{i}) ' ( j 

the last expression should be asymptotically equal to in. The results about wavepacket propagation 



give full information about the previous physical expectation values. In fact, (13) can be viewed 
as the contraction between the evolved state and a semiclassical boundary intertwiner with one 
insertion, so 



E^dim^iMii) " T, h WTMh) = *°' [ j 



where we used (12); what we have found is that dynamical and kinematical mean (asymptotically) 
coincide. We stress that, not only the peakiness of the evolved state is required in order to have the 
right expectation value, but also the phase of the evolved state must cancel exactly the phase of the 
initial intertwiner. The same properties (peakiness and right phase) hold for the 3-to-2-propagation 
(see numerical results in the next section), and the expectation value of two intertwiners turns out 
to be the correct one, i.e. in. 
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5 Numerical analysis 



We wrote a numerical algorithm performing the 4-to-l and 3-to-2 evolutions, and calculating the 



physical expectation values (13) (14); the algorithm computes very big sums serially with a method 
similar to the one in [23j|24 . The results are shown in the figures. In Fig. [3] the result of the 
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Figure 3: On the left: modulus of the evolved state for the 4-to-l propagation performed by the 
flipped vertex (jo = 30). On the right: its real and imaginary (dashed) part. 

4-to-l evolution for jo = 30 is reported. From the plot on the left (the modulus) we can see that 
the evolved state is a Gaussian peaked on i$ with the same width of the "incoming" Gaussians. 
On the right the real and imaginary parts are plotted, and it is clearly visible that the frequency 
of oscillation is — 7r/2, which is exactly the phase opposite to the one of initial packets. 
In Fig. [7] are shown the results of the 3-to-2 propagation (moduli), from jo = 10 to jo = 32 for 
even jo's. Compared with the 4-to-l case, here the Gaussian shape seems not to be conserved, 
but the state is nevertheless peaked on iq and presents a —n/2 phase in both variables; actually a 
convergence to an elliptic Gaussian is taking place (we explored up to jo = 56). Non-Gaussianity 
has to be imputed to quantum effects. Small deviations from Gaussianity are present also in the 
4-to-l-evolution, though less pronounced. Both in the 4-to-l and 3-to-2 evolution, non-Gaussianity 
gives rise to deviations of physical expectation values from the classical behavior, well visible in the 
plots in Fig. |4j 

The physical expectation values (Fig. [4]) are in complete agreement with the expected ones. The 
small deviations from the semiclassical values gradually disappear as jo increases. 



6 Correlation function 

Here we present some very preliminary results about the graviton propagator in the EPR spinfoam 
model, in the rough approximation of fixed spins. In other words, we consider the 2-point function 
over the boundary state ([6]). With our approximation, the spin-spin and spin-intertwiner corre- 
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Figure 4: On the left: physical expectation value of i\. On the right: physical expectation value of 
i\i2- The solid line is the expected behavior. 



lation functions vanish, but some components of the graviton propagator are proportional to the 
intertwiner-intertwiner correlation function and we will study these ones. The 2-point function, or 
propagator, over the boundary of a 4-simplex is defined by: 

bcd W({jh WXgg • gg - nff • ng)(gg • - n£> • n^({j}, W) 

E {jh{ii w({j},{i}M{jh{i}) { ' 

where m, n and a, b, c, d run over {1, . . . , 5}, E^ is the electric field (densitized triad) operator at 
the node m, projected along the normal n!$ in m to the face shared by the tetrahedra m and a. 
In the diagonal-diagonal components (a = b, c = d) the electric fields act as area operators, so 
that the 2-point function is essentially an expectation value of two spin insertions u Sj5j" , which 
in our case of fixed spins vanishes. In the diagonal-nondiagonal components (a = 6, c ^ d) the 
first couple of electric fields gives a spin insertion "5j" , while the second couple acts nontrivially 
(in fact the nondiagonal action is the one that "reads" the intertwiner quantum numbers at nodes) 
but also those components vanish because of the presence of the spin insertion. The only surviving 
components are the nondiagonal-nondiagonal; they are quite complicated but some of them are 
simpler, and we will consider only them. Consider in the boundary state a node m labeled by the 
virtual spin i m in a certain pairing, and concentrate only on those (ft, 6) which are coupled to im- 
As an example, if we take the node m = 1 and the surrounding spins are labeled as in Fig. |5j then 
we consider only a,b £ {2,3} or {4,5} (a ^ b). Then consider another node n labeled by i n and 
indices Cy^d coupled to i n . For those components the action of graviton operators is diagonal and 
gives insertions of the kind "(^^ — 8j — ^)(^^ — &3 ~ > so 



rabcd Z {j} , {i} W({j},{i})^n*({j},{i}) = (5 im 5i n ) 
X j 2 E {jh{i} W({j},{i})n{j},{i}) " Jo 2 ' 
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Figure 5: node 
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Figure 6: On the left: real part of the intertwiner correlation function, divided by jg. On the right: 
its imaginary part. 



If the propagator has the Newtonian scaling in the semiclassical regime, it should scale asymptoti- 
cally as the inverse of jo; equivalently, the quantity (5i m 5i n ) should scale linearly with jo. The plot 
of (di m 5i n ) from jo = 2 to jo = 50 (with step 2) is shown in Fig. |6j The scaling is clearly not the 
Newtonian one, and this could be due to our choice of boundary state, which freezes the spins, or 
maybe to some pathology of the model. The auspicious results about the evolution of wave packets 
and the physical expectation values seem to exclude the latter possibility. 



7 Conclusions and outlook 

We studied the propagation of semiclassical intertwiners over a 4-simplex, using the EPR spinfoam 
model (flipped vertex). This approach, introduced in [20], which is a way to study the semiclas- 
sical limit of spinfoam models for quantum gravity, turned out to be viable both analytically and 
numerically, and gave encouraging answers. In particular, certain coherent states turned out to 
evolve in accordance with classical general relativity. Then we read the results as physical ex- 
pectation values of observables. In the end, we showed a numerical calculation of the intertwiner 
correlation function, but the scaling law w.r.t. distance is not the one giving rise to Newton law 
in the semiclassical regime. This has to be expected, as the approximation of freezing the spins 



9 



could also prevent intertwiner fluctuations (remember that a classical 4-simplex is fully determined 
by the ten edge lengths, so if those lengths are given then the dihedral angles between triangles 
are automatically determined). Though positive, we regard our results as partial and tentative: 
one should get rid of the "fixed spin" approximation and see if the wavepacket propagation is still 
correct and then compute the 2-point function in the semiclassical limit and see if the scaling with 
distance is the Newtonian one. Further numerical [23] and analytical investigations have already 
started and we expect in the next future to be able to give more precise answers. 
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